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This paper presents an initial step toward model identification from wind tunnel data for 
an airliner configuration. Two approaches to modeling a transport configuration are 
considered and applied to both steady and large-amplitude forced-oscillation wind tunnel 
data taken over a wide range of angles of attack. Only limited conclusions could be drawn 
from this initial data set. Although model estimated time histories of normal force and 
pitching moment agree reasonably well with the corresponding measured values, model 
damping parameters did not, for some cases, have values consistent with small amplitude 
oscillatory data. In addition, large parameter standard errors implied poor information 
content for model structure determination and parameter estimation. Further investigation 
of the modeling problem for more general aerodynamic models is recommended with close 
attention to experiment design for obtaining parameters with high accuracy. 
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Nomenclature 

parameters in deficiency function 

non-dimensional aerodynamic coefficient where a = (N or m) 

non-dimensional pitching moment coefficient 

non-dimensional normal force coefficient 

mean aerodynamic chord, m 

frequency, Hz 

deficiency function 

position of aerodynamic center, % c 

nondimensional frequency 

tail arm, m 

pitch rate, rad/sec 

reference area, m 1 2 

time, sec 

velocity, m/s 

wing center-of-pressure position, % c 
angle of attack, rad 
dynamic angle of attack, rad 
initial (average) angle of attack, rad 
amplitude of angle of attack, rad 
downwash angle, rad 
tail setting angle, rad 
integration variable, sec 

time constants for normal force, center of pressure, and downwash angle 
measurement noise 
angular frequency, rad/sec 
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subscripts 


Aerodynamic Derivatives 

II II II II 
^ tel s 

amplitude 

measurement 

tail 

wing angle of attack 

= where a = N or m 

“ da 

c dC n 

C„ = where a=iVorm 

2V dq 

superscripts 


abbreviations 

B,T,W = 

body, tail, wing components 

a.c. = aerodynamic center 

BT, BW = 

body-tail, body-wing components 

d.o.f = degrees of freedom 

" = 

estimate 

c.g. = center of gravity 


I. Introduction 

R ECENTLY there has been increased interest in identification of aircraft aerodynamic models for transport 
configurations that include unsteady aerodynamic terms. One problem identified under the Aviation Safety 
Program of NASA is that of developing mathematical models that allow prediction of aircraft response over a large 
flight envelope including upset conditions. This leads to the problem of developing adequate mathematical models, 
and methods for their identification, in flight regimes with nonlinear and unsteady aerodynamic behaviors. The 
majority of models identified to date have been developed for the longitudinal motion of tailless aircraft, (see Refs. 
1-3), where the aerodynamic models had relatively simple form. In formulating aerodynamic model equations for a 
wing-tail combination, however, more complicated model forms can be expected. This problem was first addressed 
in the early twenties by Cowley and Glauert 4 . They realized that there was a time lag before an aerodynamic 
disturbance over the wing reaches the tail. They assumed that the downwash associated with a change in lift is equal 
to the corresponding steady value but the effect at the tail is delayed by the time for the airplane to travel a distance 
equal to the tail arm. The change of the lift on the tail due to downwash was included in the damping-in-pitch 
derivative. Now this change is interpreted as a tail contribution to the acceleration derivatives, i.e., the derivative of 
the lift and pitching moment with respect to rate of change in the angle of attack. 

The investigation of the downwash angle was extended by R. T. Jones and Fehlner 5 by considering both the 
growth of wing circulation and the delay in the development of lift by the tail. They presented an expression for the 
downwash indicial function associated with a change of the lift on the wing. Considerable effort to the rigorous 
explanation of the downwash both qualitatively and analytically was given by Tobak 6 . He developed indicial 
functions for the lift and pitching moment of a wing-tail combination and presented numerical results for several 
representative cases in supersonic flight regimes. The theoretical results were then applied to the longitudinal 
stability analysis of an aircraft. 

The linear aerodynamic equations for a planar motion of an aircraft with horizontal tail were developed by 
Klein 7 . Unsteady effects in these equations were expressed by linear indicial functions for wing and tail 
aerodynamics, and for the downwash angle. The emphasis was given on model structure which would explain the 
aerodynamics of wing-tail combination with sufficient accuracy and, at the same time, should be simple enough to 
provide good conditions for parameter identifiability from experimental data. 

A different model of the wing-tail combination was proposed by Khrabrov 8 et al, from the Central Aero- 
Hydrodynamic Institute (TsAGI). They assume that the aircraft model can be represented by body, wing, and tail 
components, and that the aerodynamics are, in general, nonlinear and unsteady. The resulting model included 
algebraic and differential equations with parameters dependent upon the angle of attack. The model was used in the 
analysis of steady and unsteady wind tunnel data and the identified model was presented in several graphs 
comparing measured and estimated data. Unfortunately, no numerical values of estimated parameters and their 
accuracies are given. 

The purpose of this paper is to use wind tunnel data from Ref. 8 and aerodynamic model equations based on 
those of Ref. 7 and 8 for model identification, i.e., model structure determination and parameter estimation. The 
paper is a product of ongoing research into modeling of an airliner with two lifting surfaces, wing and tail. Results 
presented in this paper are limited and should be considered preliminary. 
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II. Measured Data 


The measured data were obtained from static and dynamic testing of an airliner in the TsAGI low subsonic wind 
tunnel. The steady data included the normal force and the pitching moment of four different configurations, i.e., 
body alone (B), body-tail (BT), body-wing (BW) and the complete model (BWT). These measurements cover the 
angles of attack from -10° up to 30°. 

The dynamic data were obtained from a forced oscillatory motion of the same model in the same tunnel using the 
same test rig. The experiment was executed at different initial values of angle of attack, do, frequencies, and 
amplitudes, cxa- For aircraft model identification, large amplitude (cxa > 5°) time histories of C v and C„, for three 
configurations BT, BW, and BWT, at frequencies approximately equal to 0.4Hz, 0.8Hz, and 1.3Hz, were used. For 
large oscillation amplitudes of 15°, a 0 was equal to 5° and 15° and for oscillation amplitudes of 10°, a 0 was equal to 
5°, 10°, and 20°. Only one cycle of oscillation was given in each time history. 

The small amplitude oscillatory data were obtained for a A = 3°, a 0 from -10° to 30°, frequencies of 0.5 FIz, 1.0 
Flz, and 1.5 FIz, and two configurations, BW and BWT. The results are available in the form of the in-phase and out- 
of-phase components C Na , C Nq , C ma , and C mq (see Ref. 1 or 2) and were used for a comparison with identified 
models. The normal force components are defined as 

C Na =C Na -k 2 C Nq (1) 

C Nq =C Nq +C m (2) 

and a similar expression can be written for C ma , C mq . Because of the limited number of frequencies, the data 

cannot be used for parameter estimation of C Na , C ma , C Nq , and C mq . A more detailed description of the test facility, 
test conditions, and results are given in Ref. 8. 

III. Analysis of Steady Data 

Modeling the steady data of Ref. 8 took into account the form of measured data. It is assumed that the 
aerodynamics of the complete aircraft can be formulated as the sum of three contributions due to the body, the tail, 
and the wing. Then the corresponding models have the form 

C N {t) = C»{a) + Cl(a) + S t C T N (a t ) (3) 

C m (0 = Cl (a) + Cl Q + Cl (a)x w (a) - 1,S,C T N {a , ) (4) 


where 


a t = a + (/>-s{a) 


(5) 


S, = S, / S , £ t =i t /c , and x n is the distance between the wing center of pressure and aircraft c.g. as a fraction of 
mean aerodynamic chord. Eqns. (3) and (4) are taken from Ref. 8 with slightly different notation. 
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Figure 1. Normal-force coefficient for four Figure 2. Pitching-moment coefficient for four 
configurations (reproduced from Ref. 8). configurations (reproduced from Ref. 8). 

The values of C v and C,„ for a complete model and its components are reported from Ref. 8 in Fig. 1 and 2. in 
these plots B, BW, BT, and BWT indicate the values for the body, body-wing, body-tail, and the complete model. 

The values of Cy and C 7 ^ in Eqs. (3) and (4) were obtained from the relations 


r W _ r BW r B 
- L iV l N 


(6) 


r T _ r BT r B 


(7) 


Furthermore in Eqns. (3-5), Cjf /Q is equal to the wing pitching moment when cjj' is equal to zero. (f> is the tail 
plane setting and s is the downwash angle at the tail. From the geometry of the model S t = 0.26, l t = 4.62 and (j)= - 
10° and from C% the value of = -0.13. The variable x W (a) can be computed by removing terms associated 
with the tail in Eq. (4) and using the corresponding BW measurements. The downwash angle follows from solving 
Eqs. (3-5) for a, at selected a. The two variables x W (a) and 8(a) are plotted in Figs. 3 and 4 and compared with 

results of Ref. 8. The difference between the two pairs is noticeable for a < 5°. The reason for these differences has 
not been found. 







x w Ref. 8 
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-10 
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a 0 , deg 

Figure 3. Distance x w estimated from steady wind tunnel data. 
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15 



20 30 

a 0 , deg 

Figure 4. Downwash angle at the tail estimated from steady wind tunnel data. 


IV. Model Postulation 

Two model forms will be used in model identification from oscillatory wind tunnel data. The first model, taken 
from Ref. 8 after minor changes in the notation, will be called the TsAGI model. It has been developed under the 
following assumptions: 

a) For the one d.o.f. harmonic motion in pitch, C N = C N (a,q) and C m = C m {a,q ) , where q = a . 

b) Aircraft aerodynamics characteristics are obtained by contributions of the body, tail, and the wing. 

c) The unsteady effects in the aerodynamic model equations are included in Cj^ , x W , and s. 

Then the model equations can be formed as 

C v (t; a) = C,y (oo; a) + (t; a) + S t C^ (t\ a t ) + C Nq (go; a)d(t) (8) 

2Vq 

C m (t;a) = Cf (oo; a) + C% () + C% ( f, a)x w (f, a) - 1,S,C T N ( f, a, ) + C mq (oo; a)d(t) (9) 


where 


a t - a + </> + a d -e(t;a) 
t.a 

«.= — 

t n (a)C%s + C^s = C^ s {co;a) 
T x {a)x +x = x (oo ;a) 
T E (a)s + £ = s(w,a) 

The normal force coefficient due to the wing is partitioned as 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 
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( 15 ) 


r w _ r w r w 

- ^NA + '-'MS 

where C^ A is the linear part (attached flow) and C^ s is the nonlinear part (separated flow). Symbol oo indicates a 
variable evaluated in steady flow. Unknown parameters in the TsAGI model are the time constants, 
t n , t x , t £ , and two damping terms, C N (oo;a) and C m? (oo; a) . In general, each of these parameters can be a 

function of a. 

The second model, referred to as the LaRC model, was taken from Ref. 7 as 

C a (t) = c a (0) + \' 0 C aa (t -T-a(r))d(T)dr + ^\' 0 C ai/ (t -T-a(T))q(r)dr (16) 


where a = N or m, C a (0) is the value of the coefficient during steady conditions, and C aa (t) and C a ^ (t) are the 

indicial functions describing the response of C a (t) to a sudden change in a and q. 

Sometimes it is more convenient to use the deficiency functions, F a ^ (t) , F a (t ) , rather than the indicial 

flinctions. Then the two functions are related as 

F a a (t) = C aa (co;a)-C aa (t) (17) 


F Uq (t) = C Uq (oo; a) - C Uq (t) 


(18) 


where a = N or m. Substituting Eqs. (17) and (18) into Eq. (16) and replacing F (t) by its steady value results in 


C a if) = c a (°°; «) + ^7 C aq (oo; a)q - }' F a<x (t - r; a(r))d(r)dT ( 1 9) 


The indicial function C N ( t ) , and consequently the deficiency function F N a ( t ) , includes combined responses 

of the wing and the tail, and interference effects between those two lifting surfaces. It is therefore assumed that the 
resulting indicial functions are given as a sum of these four components: 


(1) response of C ^ to a unit step in a K , while a, = 0. 

T 

(2) response of C N to a unit step in a ln while a, = 0. 

(3) response of C ^ to a unit step in a h while a w = 0. 

T 

(4) response of C N to a unit step in a„ while a w = 0. 


For this study, however, it is assumed that the tail does not exhibit an unsteady response and that it does not affect 
wing response, therefore the third and fourth components (bottom row) are neglected. The first component 
represents the response of an isolated wing, and the second expresses the lift on the tail due to a change in the 
downwash induced by the lift of the wing. 

As pointed out in Ref. 7, the contribution of the wing and the tail to their indicial functions can be written as 

C N (t) = C^(t) + S t Cl T (a t ) (20) 

C m (t) = C mac +hC*(t)-J t S t C H J(t) (21) 

Where C,„ ic is the pitching moment about the aerodynamic center of the wing and h is the non-dimensional 

distance between the aerodynamic center and the aircraft center of gravity. For the LaRC model, Cy is located at 
the aerodynamic center. Combining Eqs. (19-21) yields 
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(22) 


c n (0 = c n (°°; «) + Cjv 9 (co',a)q(t) 
lv o 

- Io[ f aL w - *■; «(r)) + “ r; «( r ))]«( z ') t/ r 

C m (0 = c m (°o; a) + C mq (*>; a)?(t) 

2K o (23) 

“ \o[ hF N aw ^ ~ G «0")) - ~ l t^i F N T aw it ~ r;a(r))]a(r)rfr 

For model identification the deficiency functions in Eqs. (22-23) must be specified. Experience with identification of 
tailless aircraft indicates that a simple exponential function 

= ae~ bt (24) 

- /V «W V 7 

D] 

can be a good approximation. Formulation of the deficiency function F N , however, is more complicated. A 
model for the corresponding indicial function, (i ) , is presented in Ref. 7 for linear aerodynamics combining 
three effects: (1) change of C ^ followed by a change in a w , (2) change in s at the tail following a change in cj ^ , 

DT 

and (3) a change in C N due to a change in s. The resulting expression for the indicial function is 

C B J aw it) = -C B N T g (t)e(O) - J' C*l (t - T)s(T)dr (25) 

The function (t) represents the lift on the tail during passage of the tail through a step change in downwash that 
is equivalent to the passage through a sharp-edge gust 5 . 

The model for the indicial function, ( t ) , is too complicated for model identification. It can be substantially 

simplified by using the Cowley-Glauert explanation of downwash delay. According to their assumptions the 
downwash from the wing is delayed from reaching the tail by the time A t = i t IV . With this assumption, the 
indicial function can be simplified as 


c Zjt) = cB N T aw i°o)it-A t ) (26) 

V. Model Identification 

In preparation for application of the identification algorithms to measured data, all identification algorithms were 
validated against simulated data with various signal-noise levels, numbers of cycles of harmonic data, initial 
parameter error, and number of different frequencies. This analysis indicated that single-cycle harmonic data for 
each of three frequencies and modest measurement noise levels may contain limited information for identification 
and prevent convergence of the output error algorithm. In this paper the TsAGI model, given by Eqs. (8-15), was 
used in model identification. The LaRC model, represented by Eqs. (22-23), was utilized only in its simplified form 
for the BT configuration. The LaRC model has less pre -defined model structure and therefore requires more 
information content for identification. 

A. BT Model Identification 

The BT model equations were obtained by a simplification of the general equations resulting in 

C N it) = C B N {a) + S t C T N ia t ) + C T Nq ia)^ (27) 
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(28) 


CJO = C B (a)~ e t S t C T N (a t ) + C T mq (a )^~ 


for the T sAGI model and 


C N {t) = S t C BT {a t )+C BT q {a)— (29) 

C m (0 = - l t S,C BT {a , ) + C BT q ( a ) ^ (30) 

for the LaRC model. For both sets of equations 

a t = a + </> + a d . (31) 


BT BT 

Eqns. (27-30) were further simplified by assuming that the damping parameters, C Nq and C mq , do not change with 

a. After substituting measured values into the model equations, the least squares method in time domain was used 
for parameter estimation. 

The most consistent results were obtained from Eqs. (29-30) after replacing the term - H t S t C BT (a f ) by the 

directly measured pitching moment C BT (a t ) . The resulting parameter estimates are summarized in Table I which 

includes their mean values, 9 , and standard errors, er . 
The mean values shown were obtained by averaging 
estimates over 3 different frequencies for given a.„ and 

a A . The last term in Table I shows values of C BB 
computed from C m ^ = ~^t^Nq • These values are in 

agreement with direct estimates of C m ^ . All parameters 

were estimated with large standard errors, i.e., low 
accuracy, suggesting limited information content in the 
data. 

Estimated parameters are also compared with 
derivatives C Nq and C mq in Figs. 5 and 6. These derivatives were estimated from small amplitude oscillatory data 

at the frequency of 1.5FIz. For a comparison of different parameters it is assumed that the difference between C Nq 

DT D T 1 

and C Nq , and C mq and C„^ are caused mainly by parameters C N& and C m ^ or their nonsteady counterparts (see 
Eq. (19) and Ref. 9). The direct comparison of parameters in both figures is not possible but at least some trends in 

DT 

the results can be observed. The parameters in Fig. 5 indicate compatibility between C Nq and anticipated values of 
C Nq using the relationship C Nq = C Nq - C N & . A similar argument can be made for the differences between C m ^ 
and C mq in Fig. 6. 
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Table I. Estimated damping parameters from 
large-amplitude oscillatory data, BT configuration. 


an, deg 

a A , deg 

r BT 
L Nq 

^BT 
^ mq 

-4 

r BT 
L Nq 



9 

(O') 

9 

(a) 

9 

(a) 

5 

15 

3.1 

(0.82) 

-15 

(1.7) 

-15 

(3.3) 

15 

15 

2.9 

(0.99) 

-15 

(2.1) 

-13 

(4.6) 

5 

10 

3.0 

(0.98) 

-12 

(2.5) 

-14 

(4.5) 

10 

10 

2.7 

(1.40) 

-12 

(1.4) 

-12 

(6.5) 

20 

10 

1.5 

(0.49) 

-17 

(2.6) 

- 7 

(2.3) 
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Figure 5. Estimated normal-force damping 
coefficients due to pitch rate from oscillatory 
data for BT and BWT configurations. 



Figure 6. Estimated pitch damping coefficients 
from oscillatory data for BT and BWT 
configurations. 


B. BW Model Identification 

The BW model equations used in data analysis were taken from the TsAGI model. These equations are 


C N (t) = C b n (oo; a) + Cl (a , ) + C«f («; a) ^ 


(32) 


C m (0 = Cl (oo; a) + Cl 0 + Cl (a)x w {a)+ C™ (oo; a) ^ 


(33) 


As in the complete model, the normal force coefficient was partitioned as 


r w _ r w , r w 
c JV “ L NA + L NS 

then a time delay due to flow separation can be estimated from 


(34) 


t n (a)Cls + Cl s - Cl s (oo; a) 


(35) 


Similarly, for identification of the pitching-moment model, the coordinate of the wing center of pressure can be 
estimated from 


t x (a)x w + x w = x w (oo; a) (36) 

Now the estimation problem for the C N model parameters can be formulated as follows: given time histories of 
C N (t) , a(t) , and a(t) , and analytical functions Cl s (x;a), Cl A (co;a) , and cl(co;a), estimate unknown 

constants t n and . To simplify notation, introduce r/ = cl s , 0\ = r, v , 0 2 =-^-C^, u = Cl s (co;a) , 

y = C N it) -Cl(co;a)~ d A (co;a), then 
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6 x r] + 77 = u 


(37) 


y = rj + 0 2 a 


(38) 


A similar formulation can be made for C m where the unknown parameters to be estimated are r x and C^j . The 

unknown parameters in the state-space models were estimated by an output error method explained in Ref. 10. 

Parameter estimates and their 
standard errors are given in Table II. 
Both time constant and damping 
parameters were considered independent 
of a for each data set analyzed. The 
average value of r N ~ 13 for all five 
cases agrees with limited previous 
experience. The estimates of the second 
time constant exhibit large scatter around 
the average value of r x =41. This value 
could not be verified. Also the effect of 
z x on estimated C m (t) was not 
investigated. Large variability of mean values and large standard errors indicates the low accuracy of parameter 
estimates and limited information content of the measured data. 

Damping parameters in Table II were compared in Figs. 7 and 8 with the estimates from small amplitude 
oscillatory data. Using the same assumptions discussed with Figs. 5 and 6 for comparing the damping estimates with 
out-of-phase measurements, an agreement between C Nq parameters can be argued for a from -5° to 5°, and for a 

near 25°. A large discrepancy in C Nq near a=15° has not been explained. Small amplitude oscillatory data in Fig. 8 
suggests very small negative values would be appropriate for C mq and confirms that the positive values estimated 
for C mq are incorrect. 


Table II. Estimated damping parameters and time constants 
from large-amplitude oscillatory data, BW configuration. 


a 0 , deg 

a A , deg 

t n 

r BW 
L Nq 

Tt 

r BW 



e (<x) 

e (a) 

e (&) 

e (<r) 

5 

15 

12.5 (1.2) 

4.8 (2.2) 

50.0 (3.2) 

5.1 (.63) 

15 

15 

9.6 (1.4) 

8.8 (3.7) 

40.6 (2.6) 

7.4 (.78) 

5 

10 

14.4 (1.9) 

5.8 (2.9) 

25.6 (1.7) 

3.8 (.71) 

10 

10 

12.7 (1.6) 

10.5 (4.0) 

34.0 (1.6) 

10.5 (.69) 

20 

10 

16.6 (1.4) 

13.7 (4.4) 

56.2 (.67) 

0.2 (.67) 




Figure 7. Estimated normal-force damping Figure 8. Estimated pitch damping coefficients 

coefficients due to pitch rate from oscillatory from oscillatory data for BW configuration, 

data for BW configuration. 

Measured and computed time histories of C N and C m using parameter values from Table II are plotted in Figs. 9 
and 10. Although forced-oscillation experiments were performed separately for each frequency, time history plots 
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show all 3 frequency cases stacked into one plot for analysis. In general, there is good agreement between these 
plots. Some effect of uncorrected model error in C m is visible in part of the data with f = 0.4Hz. Figs. 11 and 12 
show both coefficients plotted against a with f = 0.8 Hz. Included in both figures are also steady values of both 
coefficients. Part of the C N (a) curve for 10° < a < 20° with pronounced unsteady effect is captured quite well by 
identified model despite inaccuracies in both damping parameters. Unexplained differences, however, are apparent 
in the linear part of C N (a). The correlation between measured and computed C m (a) is poor. Finding an 
explanation for this requires additional study. 



time, sec 


Figure 9. Measured and estimated time histories 
of normal force for BW configuration, a»=5°, 
and a A = 15°. 



time, sec 


Figure 10. Measured and estimated time 
histories of pitching moment for BW 
configuration, a 0 =5°, and a A = 15°. 



Figure 11. Variation of measured and estimated 
normal force with a for BW configuration, 
a 0 = 5°, and a A = 15°. 



cc 0 , deg 

Figure 12. Variation of measured and estimated 
pitching moment with a for BW configuration, 
a 0 = 5°, and a A = 15°. 
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C. BWT Model Identification 

For BWT model identification only the TsAGI model was considered. It is represented by Eqs. (8-15). The 
unknown parameters in these equations are the two damping derivatives and the time constant, r £ , characterizing 
time dependency of downwash at the tail. The remaining terms in the model equations were obtained from steady 
data and estimates from BT and BW model identification. As in the previous cases, it was assumed that for each 
data set the unknown parameters are independent of a. 

Unknown parameters were estimated using three cycles of oscillatory data combined into one set for a given a 0 
and a_ 4 . Each cycle in the set was measured at a different oscillation frequency. Estimates of t £ were consistent with 
an average value of 26. The parameter estimates for C Nq varied from -5 to 4 and for C mq from 14 to 24. The 
negative values for C Nq and positive values for C mq can’t be substantiated. They are also in disagreement with 

estimates from BT model identification as seen in Table II. 

Figures 13 and 14 compare measured and computed coefficients C N (t) and C m (t ) . The measured data were 
taken from three cycles at a 0 = 20° and a A = 10°. The computed coefficients were based on parameter estimates 
t £ = 26, C Nq = -5, and C mq = 14. The agreement between measured and computed coefficients is very good. To 
investigate the contribution of the damping terms in this case, time histories of C N and C m were recomputed for 
the same magnitude of damping parameters but with opposite sign, i.e., C Nq = 5 and C Nq = -14. The resulting 
values of C N (t) and C m (t) were almost the same as in the preceding computation. These results indicate very low 
sensitivity of both coefficients, C N and C m , to damping parameters for this model and data set. At present, this 
problem has not been explained. It is expected that further investigation of this problem will continue together with 
experimental design for obtaining parameters with high accuracy. 




Figure 13. Measured and estimated time Figure 14. Measured and estimated time 

histories of normal force for BWT configuration, histories of pitching moment for BWT 

a 0 = 20°, and a A = 10°. configuration, a 0 = 20°, and a A = 10°. 


VI. Concluding Remarks 

This paper presents aircraft model identification from wind tunnel data of an airliner and its components: body, 
body-tail, body-wing, and the complete model. The measurement data included results from static and oscillatory 
tests over a wide range of angle of attack. Two mathematical models of an aircraft with one degree-of-freedom in 
pitch were postulated. These models differ mainly in the formulation of unsteady aerodynamics. 
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The measured steady data were fitted by high order polynomials in angle of attack. From these polynomials 
contributions of the tail and wing alone to normal force and pitching moment were computed followed by 
determination of downwash angle at the tail and center of pressure position of the wing. 

Parameter estimation started with the analysis of the body-tail oscillatory data. After substituting measured data 
into the model equations for the normal force and pitching moment coefficient, damping parameters in these 
equations were estimated by a least squares method. The model for the body-wing configuration included four 
unknowns: two damping parameters and two time constants. These time constants characterize the dynamics defined 
by differential equations for the unsteady normal force and center of pressure position. After reformulating model 
equations into a state-space form, the parameters were estimated by an output error method. In the complete model 
the two time constants from the body-wing analysis were treated as known parameters. The number of unknowns for 
the complete model analysis was thus reduced to two damping parameters and one time constant for the downwash 
angle dynamics. 

All models with estimated parameters fit the measured data quite well. The estimated time constants for the 
normal force and downwash were consistent with small deviations from their average values. On the other hand, 
estimates of the time constant for center of pressure position exhibited large scatter. The value of all three time 
constants couldn’t be verified because there are no references dealing with the same estimation problem. 

The estimates of damping parameters of the body-tail and body-wing configurations were compared with the 
out-of-phase components estimated from small amplitude oscillatory data. The damping parameters agreed with a 
trend in the out-of-phase results with the exception of the wing contribution to the damping in pitch parameter where 
those values were positive. 

The estimates of both damping parameters in the complete model were mostly unexpected values with positive 
signs. This is, of course, unacceptable for an aircraft with wing-tail combination operating in the post-stall 
conditions. It was demonstrated by simulation that both normal force and pitching moment models were insensitive 
to both damping terms and this was reflected in the identification problem for this test data. Insensitivity of model 
parameters and large standard errors suggest low information content in the data. It is expected that the ongoing 
research will explain the problems that occurred in this study and also address the experiment design problem for 
obtaining parameters with high accuracy in more general aerodynamic models. 
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